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Abstract. Charmonium formation via charm quark in-medium recombination in heavy ion inter- 
actions at collider energies has the potential to probe some properties of the medium by utilizing 
the sensitivity of the recombination process to the momentum distribution of the quarks. We have 
examined the transverse momentum spectra of J J y/, characterized by (pr 2 ), which result from the 
formation process in which the charm quark distributions are unchanged from their initial produc- 
tion in a pQCD process. This is contrasted with the case in which the charm quarks have completely 
come into thermal equilibrium with an expanding medium whose properties are determined by the 
spectra of produced light hadrons. We find that the resulting (pr 2 ) of the formed J / y/ provide a dis- 
tinct signature of the underlying charm quark spectra, and that signature is essentially independent 
of the detailed dynamics of the in-medium formation reaction. In addition, both of these signatures 
are sufficiently separated from the case in which no in-medium formation takes place. Finally, uti- 
lizing a model for the fraction of J/\]f which originate from in-medium formation, we predict the 
centrality behavior of these signatures. 



The role of J / produced in high energy heavy ion collisions as a signature for color de- 
confinement [ 1J] has evolved in recent years with the realization that at collider energies 
an additional formation mechanism may become significant M,\M- This depends on the 
initial production of multiple cc pairs in sufficient numbers. Initial estimates [0] from 
extrapolation of fixed-target data put this number at about 10 for central Au-Au colli- 
sions at RHIC. Subsequently, measurements by the PHENIX and STAR experiments 
indicate even higher numbers, between 20 [5] and 40 |6], respectively. The in-medium 
formation picture we consider here uses competing formation and dissociation reactions 
in a Boltzmann equation to calculate the final J/\y population. The absolute value of 
this formation was found to be very sensitive to the underlying charm quark momen- 
tum distributions CD. In addition, there is significant dependence on details of the size 
and expansion profile of the deconfinement region, for which various model parameters 
must be introduced. The initial PHENIX data [8] suffered from low statistics, and was 
compatible with a fairly large region of model parameter space B9D. 

Recent work in this area concentrated on finding a signature for in-medium J /\y for- 
mation which is independent of the detailed dynamics and magnitude of the formation. 
We found that the pj spectrum of the formed J/\jf may provide such a signature iflflll . 
Our first calculations of in-medium formation used initial charm quark momentum dis- 
tributions from NLO pQCD amplitudes to generate a sample of cc pairs, supplemented 
by an initial-state transverse momentum kick to simulate confinement and nuclear ef- 



fects. In the absence of in-medium formation, pjof the J/\y follows from that of the 
initially-produced individual cc pairs, but one must use a model for the magnitude of 
the hadronization process. For the normalized pr spectrum, we start with that of the pair 
spectrum. In the evolution of the interacting system size from pp to pA to AA collisions, 
the prwill be in general increased due to initial-state effects of interaction of constituents 
in the nuclei. One can express this effect as 

(PT 2 ) P A-(PT 2 ) P P =A 2 [n A -l], (1) 

where is the impact-averaged number of inelastic interactions of the projectile in 
nucleus A, and A 2 is the square of the transverse momentum transfer per collision. For 
a nucleus-nucleus collision, the corresponding relation is 

(Pt 2 )ab-(pt 2 ) pp = A 2 [n A +n B -2]. (2) 

The PHENIX measurements of J/y pyspectra in pp and minimum-bias d-Au interac- 
tions allow us to determine the amount of initial state kj needed to supplement 
our collinear pQCD events. One can then extrapolate to Au-Au and predict the spec- 
trum of J/xjr which are produced from hadronization of the initial "diagonal" cc pairs, 
again for minimum bias interactions. (We use diagonal to distinguish these pairs from 
the "off-diagonal" combinations which contribute to in-medium formation.) One 
finds h~A = 5.4 for minimum bias d-Au interactions at RHIC energy (using o pp = 42 
mb), which leads to A 2 = 0.35 ±0.14 GeV 2 . We note that the relatively large uncer- 
tainty comes entirely from the difference in pj broadening in d-Au between positive 
and negative rapidity. 

Our prediction for the "normal" evolution of the pr spectrum in Au-Au interactions 
is shown by the triangular points in Fig[T] There is some scatter at high p^due to the 
limited statistics of the number of pQCD-generated pairs. The (pr 2 ) of this spectrum is 
approximately 6.3 GeV 2 . 

The properties of the normalized pt spectrum for the formation process follow from 
two separate effects: First, the fact that the process is dominated by the off-diagonal 
pairs introduces a modified initial /^distribution. Next, one weights these pairs by a 
formation probability for J/\j/. We use the operator-product motivated cross section 



1 13l 11411 for cc forming J/\]/ with emission of a final- state gluon, which of course is 
just the inverse of the dissociation process. However, any cross section which has the 
same general properties as this one gives essentially the same result [10]. We show by 
the square points in FigU]the prediction for the formed J / 1//. One sees that this spectrum, 
characterized by (pr 2 ) approximately 3.6 GeV 1 , is substantially narrower than the one 
with no in-medium formation. 

We next considered charm quark momentum distributions which would follow if the 
charm quark interaction with the medium were so strong that they come into thermal 
equilibrium with the expanding region of deconfinement. The parameters of temperature 
and maximum transverse expansion rapidity are determined by a fit to this thermal 
behavior of the produced light hadrons. The application to charm quarks was originally 
motivated in Ref. lfl5ll . who showed that the low -pr spectrum of decay leptons from 
charmed hadron decays would not be able to differentiate between the thermal and a 
purely pQCD distribution. We show here, however, that the pj spectrum of in-medium 
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FIGURE 1. Comparison of im-medium J/\j/ transverse momentum spectra predictions from various 
scenarios. 



formed J/y is very sensitive to this distribution flfjll . The circles in Fig|T]result from 
formation calculations using T = 170 MeV and yj ma x = 0.5 for the thermal charm 
quarks. One sees that this pj spectrum is narrower yet than in-medium formation from 
pQCD quarks, with (pr 2 ) approximately 1.3 GeV 2 . Finally, we show by the stars the 
prspectrum of J /iff which themselves obey this thermal distribution. The resulting 
spectrum falls between the in-medium formation spectra for either pQCD or thermal 
charm quark momentum distributions, with (pr 2 ) approximately 3.0 GeV 2 . 

We now proceed to investigate the variation of the pQCD-based results with respect 
to the collision centrality in Au-Au interactions. First, we use the value of A 2 extracted 
from pp and pA data, together with values of Ra calculated as a function of collision 
centrality, to recalculate the (pr 2 ) values for either the initial production or the in- 
medium formation separately. This provides the centrality behavior of the J/yf spectrum 
in the case that one or the other of these mechanisms is solely responsible for the 
total J/\jf population. We show these results together in Fig. |2j One sees as expected 
that (pr 2 ) is maximum for the most central collisions, but the absolute magnitudes 
are widely separated for initial production and in-medium formation at each centrality. 
One should note that the uncertainties are dominated by the difference between the pj- 
broadening measurements at positive and negative rapidities in the d-Au interactions. 
Thus the point-to-point uncertainties are much smaller for the centrality behavior. We 
have also included separate values for {pr 2 ) in the region limited by pj< 5 GeV, to 
facilitate comparison with experiment in this same range. 
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FIGURE 2. Centrality dependence of J/\i/(pt 2 ) contrasting predictions assuming either 100 % pro- 
duction from initial cc pairs or 100 % in-medium formation. 



In order to provide a meaningful prediction for the overall J/yr spectrum, one should 
of course include both initial production and in-medium formation together as sources. 
This requires some estimate of the relative magnitudes of these processes, and is subject 
to considerable model uncertainties. What we can say, however, is that in-medium 
formation will be most dominant for central collisions, where the quadratic dependence 
on N C c is enhanced. Conversely, one expects that initial production will increase in 
relative importance for very peripheral collisions. To get an approximate idea of how 
this effect will appear, we revert to our original model calculations which included the 
absolute magnitude results [7]. The relevant parameter is the number of initial pairs, N C c, 
parameterized by its value at zero impact parameter. These results are shown in Fig. |3] 
for two representative values of N C c =10 and 20. The anticipated dominance of each at 
opposite ends of the centrality scale is seen to be realized. In addition, since the centrality 
dependence of each separate contribution is in the same direction, the dependence of the 
total (pr ) is somewhat more flat than either of them separately. 

Finally, we include a graphical representation of the evolution of (pr 2 ) as a function 
of system size in FigurelU For simplicity of presentation, we revert back to the minimum 
bias case. In the absence of in-medium formation, the behavior of (pr } is monatonic 
increasing with system size. Including in-medium formation as a significant component 
for minimum bias collisions leads to a reduced (pr 2 ) in Au-Au interactions. The nu- 
merical results indicate that this reduction is such that the predicted Au-Au value is 
even below that measured in d-Au, i.e. a non-monatonic behavior. Hence we claim that 
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FIGURE 3. Centrality dependence of J/\I/(pt 2 ) using a model calculation to estimate the relative 
contributions of initial production and in-medium formation. 



observation of such a non-monatonic behavior can be taken as a signal of in-medium for- 
mation, while the absolute value of such a behavior will be correlated with the specific 
charm quark distribution in the medium. 
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